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ABSTRACT: Even the elusive neutrinos are trapped in matter, albeit transiently, in several 
astrophysical circumstances. Their interactions with the ambient matter not only reveal the 
properties of such exotic matter itself, but also shed light on the fundamental properties of 
the neutrinos themselves. The physical sites of interest include the early universe, supernovae, 
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and newly-born neutron stars. Detection of neutrinos from these vastly different eras using the 
new generation of neutrino detectors holds great promise for enhancing our understanding of 
neutrino-matter interactions and astrophysical phenomena. 
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1 INTRODUCTION 



In accelerator neutrino physics, a neutrino interacts once, with a constituent of 
the matter in a detector, to produce a signal. But in astrophysical situations 
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the surrounding matter can make a big difference to the neutrino physics. Some 

examples are: 

(1) In their passage from the core of the sun to the solar surface, electron neutrinos 
(f e ) feel an "index of refraction", or altered energy-momentum connection, that 
is different from that for the other neutrino species, v^iVr- This difference is 
induced by forward scattering from electrons in the solar plasma. If there is 
in addition a "mass matrix" that would cause vacuum v flavor oscillations, the 
changing index of refraction seen by a v e as it moves into regions of less density 
can give rise to the MSW effect. Q 

(2) The core of a collapsing star, and subsequently the interior of a newborn 
neutron star, or proto-neutron star (PNS), formed during the supernova pro- 
cess, is opaque to neutrinos. Initially, because neutrinos are trapped in the PNS, 
deleptonization of matter has yet to occur and the total number of leptons per 
baryon is about 0.4. Only after times of tens of seconds can neutrinos diffuse 
outwards and escape, but in so doing, they heat the matter through which they 
pass in a process reminiscent of Joule heating. Following this deleptonization is a 
cooling epoch, during which neutrinos continue to transport energy to the star's 
surface. The combination of extreme neutrino degeneracy and high temperatures 
and densities creates a unique environment in which the detailed neutrino pro- 
cesses and the resulting transport phenomena play an essential role. The escaping 
neutrino fluxes are crucial for the supernova dynamics, and are also important 
in the heavy element nucleosynthesis that might occur in the ejected envelope of 
the formerly collapsing star. Certainly, they determine the characteristics of the 

neutrino pulse that can be observed in terrestrial detectors (|2|). 
1 A recent review of neutrino oscillation physics is given in Ref. ([[]). 
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(3) In the standard model (three v flavor) evolution of the early universe, in the 

temperature region of, roughly, 0.2 < ksT/MeV < 10, the v's and P's are very 

nearly in thermal equilibrium with the other particles in the plasma until the 

temperature drops below 1 MeV or so. The corrections to the distributions near 

this point of v decoupling have been calculated quite well and provide a small 

change in the predictions of primordial 4 He abundance. But there are models, 

involving flavor mixing and the possibility of a fourth neutrino flavor, in which 

the evolution of neutrino distributions, including the effects of interactions with 

the surrounding plasma, must be studied carefully at higher temperatures (see 

Section 4 below). 

1.1 Orders of Magnitude 

We choose units so that fi = c = ks = 1. Then the weak coupling constant is 
given by Gp = 1.166 x 10~ n (MeV) -2 . Number densities in these units for some 
systems of interest, are: 

Solar center (electrons); n e w 6 x 10 25 cm~ 3 «5x 10~ 7 (MeV) 3 
Earth core (electrons); n e ~ 5 x 10 -6 (MeV) 3 

Supernova-core to neutrinosphere (nucleons) = [10 4 — 10 6 ] (MeV) 3 
Early universe in the region, 1 < T/MeV <100; n part i c i es « 5 x T 3 . 

We can estimate the length scales over which the effects of neutrino interactions 
with the media can be appreciable in the various cases: 

(a) At solar center densities, the index of refraction correction from v e -e forward 
scattering introduces a phase change of 2ir in a distance 
f« {n e G F )- 1 w 4 x 10 18 (MeV)" 1 » 2 x 10 7 cm. 
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(b) The solar center mean free path for v scattering is: 

X f « (n e G\ElY x « 2 x 10 28 (MeV)- 1 «4x 10 17 (MeV/^) 2 cm. 
This is far too long to be of any interest. 

(c) The supernova center mean free path for neutrino scattering (but subject to 
substantial medium-dependent corrections to be discussed in this review) is 

Xf « (n N G 2 F El)- 1 « 2 x 10 5 (MeV/-E„) 2 cm. 

(d) For the early universe it is the neutrino collision time that is of interest 
A//cw W- 21 /(G 2 F T 5 ) « (T/MeV)- 3 s. 

This, coupled to the expansion time, i(s)« [T(MeV)]~ 2 , shows that i/'s decouple 
from the medium at T ~ 1 MeV. 

1.2 Scope of this Article 

In this review, we shall concentrate on the physics of neutrinos in supernovae, 
PNSs and in the early universe. For these problems one needs to follow the 
evolution of neutrino distributions in time. In all cases, we envision beginning 
with an initial condition with postulated or known neutrino distributions. In the 
early universe problem, the initial time can be taken to be sufficiently early for u 
oscillations not to have affected the distributions; an initial temperature of tens 
of MeV suffices in the models that will be discussed. 

In the supernova and PNS problems, we focus on the microphysics of neutrino 
interactions with matter. For the supernova case, neutrinos are initially freely 
streaming and act mainly as a cooling mechanism. Shortly after collapse begins, 
the neutrino mean free path becomes less than the size of the collapsing core, 
and neutrinos become trapped on dynamical time scales. During this period, 
there is a complex interaction between transport and hydrodynamics. While 
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we do not discuss supernova simulations in detail, we summarize the literature 

of recent results. For the PNS problem, which is quasi-hydrostatic, the initial 

distribution of neutrinos is that of a highly degenerate Fermi gas. We present 

results of transport simulations in PNSs for a variety of assumptions regarding 

the composition of matter. 

In all three problems the dynamics of the system provides a time changing 
environment. And in these problems the particles other than v's in the soup stay 
in near thermodynamic equilibrium. 

We attempt to pose cleanly some of the questions that have occupied theorists 
in these areas, to sketch somewhat qualitatively the methods that have been used, 
and to give some examples of results. A prototype neutrino evolution calculation 
involves three steps: 

(1) deriving a kinetic equation that describes the macroscopic development of the 
evolution of the neutrino distribution; 

(2) determining rate functions in the kernel of this equation that describe the 
local scattering, production and annihilation of the neutrinos; and 

(3) solving the equation for a physical configuration. 

Calculations of phenomena in this area are necessarily approximate but we 

shall attempt to define precisely the quantities that enter the equations. We 

must caution the reader that in most of the problems that we shall discuss 

there are, at present, few definitive results. In the first place, to the extent 

that neutrino mixing (or "oscillations") are involved, the mixing parameters are 

still undetermined. Furthermore, it is essential to follow the time development 
2 In the early universe, as we cool from, say, T = 10 MeV to 0.2 MeV the neutrons and 

protons go far out of chemical equilibrium, but the nucleons are far too sparse to affect the 

neutrino distributions significantly. 
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of non-equilibrium distributions of several species simultaneously. Finally, in 

the supernova and PNS applications, nuclear forces play an essential role in the 

opacities, and the macroscopic environment to be used in the supernova case is 

dependent on hydrodynamical calculations as well. 

We warn the reader that the references are not exhaustive. 

2 NEUTRINO MIXING AND INDEX OF REFRACTION EF- 
FECT 

The basis for what follows will be the standard model of neutrino interactions, 
in which all coupling is to the left handed neutrino current operators, in places 
supplemented by neutrino mass and mixing effects, or by the inclusion of a "ster- 
ile" neutrino, which has no interactions outside of the mixing term. If there is 
neutrino mass and mixing the standard model is augmented with a term of the 
form H m i x = J2 a ,/3 m a,/3'tpaipi3, where the ip's are the v fields and (a, (3) are flavor 
indices, taking the values e, fi, r and, in the case of the consideration of a sterile 
neutrino, s. 

The differences in the forward scattering amplitude of the different flavors 
of z/s in the various media, as calculated in first order in Gp, give rise to a 
relative energy shift matrix A Q)( g that is first order in Gj?, independent of the 
neutrino energy, diagonal in the flavor indices, and proportional to the density 
of scatterers in the medium. In the usual MSW considerations this comes from 
the electron density, n e , only and is given by A aj/ g = 6 a:e 5p :e V2GFn e . However, 
in the interior of a supernova there may be enough fi~ present briefly to make 
a significant contribution to the forward scattering differences. Sterile neutrinos 
would add to the matrix, A, as would one loop weak corrections under some 
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circumstances (||). 

In all cases of interest, we have m, A <C E v . Then standard Dirac technol- 
ogy combines the neutrino mass and index of refraction effects in an effective 
Hamiltonian for the v, v system, not including interactions, except for forward 
scattering, 

h {v) = E {4(pMp) 

where E = |p| and the A matrices are given by X^\E) = m 2 /(2E) ± A. Here 
a Q (p) and b a (p) are the respective annihilation operators for left-handed neu- 
trinos and antineutrinos of flavor a. The small admixtures of right handed v's 
created by the mass term, for the case of Dirac neutrinos, are inconsequential in 
our applications. The magnitudes of the elements of the mass 2 matrix that have 
been suggested in the literature range from 10~ 15 (eV) 2 to l(eV) 2 depending on 
the application. 

3 KINETIC EQUATIONS 

If there are no v oscillations, the kinetic equations will be for distribution func- 
tions, p a (p, r )t) for each flavor a of v and P, where the r dependence is macro- 
scopic in scale and originates from spatial variation of the properties of the 
medium, or in the initial condition taken for the neutrino distribution. For the 
case with neutrino oscillations, however, there can be off-diagonal flavor coher- 
ence that extends over large distances, and the system is no longer described in 
terms of a set of probabilities for the occupancy of the modes of the neutrino 
fields for each flavor. Instead we introduce a v momentum-flavor density matrix, 
and discuss its evolution. 

At t = we take the density matrix for the complete system to be of the form of 



ES, 



+^(P)MP) E6 a> p+\ { J(E) }, (1) 
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a product of a neutrino matrix times an everything-else matrix, the everything- 

else being in thermal equilibrium; of course, these become entangled at later 

times. Introducing the operator, A(k, t = 0) aj/ g = aT(k) a a(k) J g, we define the 

momentum-flavor density matrix, 

p a AKt) = Z£ Tr[e- ( ^-^^ )/T ^^(i ,u's\A(k,t) a ^\i y 8 )] , (2) 

i 

where A(k, t) = exp(iifi)A(k, 0) exp(— iHt) , H^ u is the Hamiltonian with all 
terms containing the v coordinates left out, and Z^ v is the partition function 
associated with this Hamiltonian. The states, \i , v's), represent all possible multi- 
v states, indexed with i. The coefficients Wi in the inner sum might be chosen 
in ways that are far from a thermal equilibrium configuration (for example, such 
that there are initially no sterile i/'s, in models that contain a sterile v). H is the 
complete Hamiltonian. We take the neutrino interaction with other particles to 
be of the formP] 

H int =dJ2 Ca,p^a [ifi x [other fields] , (3) 

where the matrix £ acts in flavor space, and is generally diagonal in flavor. For 

example, if we are discussing active-sterile v mixing we have, Ca,a = 1> Cs,s = 

Ca,s = Cs,a = 0. For v e mixing with an active species, a, in an environment 

containing electrons and positrons but no other leptons, we have, £ e)e = 1 + c e , 

(a,a = 1) where c e gives the effect of the charged current e, v e scattering term. 

There is no exact equation governing the time evolution of the function, p a ^(k, t). 

Indeed, the evaluation of the two neutrino fields in A(k) at the same momentum 

k rules out the derivation of such an equation. But if we had defined a density 
3 Generalizations to a sum of terms of this form are immediate, but give lengthier expressions 

in the evolution equations. 



10 Prakash, Lattimer, Sawyer, and Volkas 

matrix with momentum indices, k, k' then, at times large compared to a typical 

EjJ 1 we would find that, for all practical purposes, this density matrix was very 

nearly diagonal in these indices. We can qualitatively describe the simplification 

as due to the rapid decoherence between the parts of the v wave function that are 

scattered at appreciably different momentum transfer (and/or energy transfer) 

from components of the medium at different (randomly situated) points in space, 

the density matrix in momentum space becoming nearly diagonal so fast that 

we need to consider only the diagonal elements in the p space. Finally, we 

consider only systems in which the local averages of the densities, compositions, 

and temperatures varies slowly in space, compared to all quantum mechanical 

distances except for the neutrino oscillation lengths. The evolution equation is 

then, 

J^(k, r, t) + |k|- x k • Vp(k, r, t) = -i[X(E),P(k, r, t)] + 
+\ E [Cp(ki, r, i)C[l - p(k, r, t)\ + [1 - p(k, r, t)Rp(ki, r, t)(] r(k x ,k) 

ki 

~\ E [Ct 1 - P( k i> r < *))Cp(k, r, t) + p(k, r, t)C[l - p(k!, r, t)](] T(k, ki) . (4) 

ki 

The function, T(k, ki), is the differential rate at which neutrinos would scat- 
ter from energy k to momentum ki if we took £ = 1, that is, if the neutrino 
interaction were flavor independent. | 

Many alternative paths to a derivation of equations like (||) are found in the 

literature. If we consider only the terms that are linear in p, as for the case 
4 As a concrete example, in the case in which one of the neutrinos is an electron neutrino, 

the plasma is constituted of electrons and positrons, and we consider v e , v^, mixing, F stands 

for rate of scattering of u^, and £ is given as above. The squared amplitude contributions from 

scattering from electrons get added to those from positrons. So we can use <^ ee = 1 + c, just as 

though there were a single term in equation (W) . 
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of a single neutrino traversing the medium, then with the right identifications 

equation (|5|) can be seen to be the same as the "quantum kinetic equation" 

of McKellar and Thomson (|j), taken in the same limits. Using a somewhat 

idealized form of interaction, Raffelt, Sigl, and Stodolsky (|5|) have given an elegant 

derivation that includes Fermi statistics (i.e. the nonlinear terms on the RHS); 

they call the result the "non-Abelian Boltzmann equation." Again for the energy 

conserving case, but without Fermi statistics, Loreti and Balantekin (|6|) have 

used the equation, which they describe as a variant of the "Redfield equation" 

to analyze MSW transitions in media with irregular electron density variations 

in space. In the vast literature from the condensed matter and quantum optics 

side, such equations are often called "Bloch equations." Typical derivations are 

given in Refs. A simple derivation of equation (|j) that captures most of 

the essential aspects of the physics, for the case of the density matrix of a single 

neutrino, is given in Ref. (||). 

A simple application of equation (Q) is to active-sterile mixing, involving two 

neutrino states that mix, only one of which scatters. For a translational and 

rotationally invariant system, we calculate a density matrix p(E, t). We keep only 

the linear terms for illustration, and we adopt the representation p = ^(Pq+P-ct), 

A = V • a . The projection operator on the active state is £ = (1 + cr z )/2, where 

ex's are the Pauli matrices. Substituting in equation (||) gives 

jW, t) = V(E) x P(E) - I [xP x (E) +yP y (E)] £ T(E, E 1 )+z^-P (E, t) (5) 



dt 
and 



8 P (E,t) = -\[P z {E)+P {E)}Y j nE,E l ) + \Y j nEuE)[P z {E l ) + P {E l )] 



(6) 
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For later use in Section (4) , it is helpful to define a decoherence function through 

D = (££ i r(E', Ei))/2 and a repopulation function through R = dPo/dt. 



3. 1 Extensions of the Evolution Equation 

In the kinetic equation (||) the sources of change of the neutrino distribution are 
transport, flavor precession, and scattering. We sketch briefly the other effects 
that should be included. For economy in notation we can consider the spatially 
uniform case, in which the v density matrix is given by P a ^(E), as in equation (|5|). 

A. vv Pair Production and Annihilation: These processes occur due to 
the conversion of charged lepton-anti-lepton pairs, or from v pair bremsstrahlung 
processes. To include these, we can introduce a set of densities for antineutri- 
nos and incorporate the coupling between the two sectors cominmg from pair 
processes. However, it is neater, to define the function [1 — p{— E)] as the D 
density matrix, to use the extension of \{E) to negative energies provided by 
the definition in equation @, and to extend the functions, V(E, E\), so that 
for E < 0, Ei > it is the pair production rate, for E > 0, E\ < the pair 
annihilation rate, and for E < 0, E\ < the anti-particle scattering rate, but 
from energy, —E\, to energy —E. With these conventions the only change to be 
made in equation (||), to accomodate pair creation and annihilation will be to 
relabel the energy sums in equation (|5|) to be Hbi=-oo- I n the case A = 0, ( = 1, 
the system would move toward a thermal equilibrium described by a v, v Fermi 
distributions with a single chemical potential determined by the excess lepton 
number for each flavor. 
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B. Expansion of the Universe: This is taken into account by using the 

time dependent temperature in the statistical factor and adding to the LHS of 
equation (||) the term — [a(t)/a(t)]E(d/dE)p(E,t) where a(t) is the scale factor. 

C. Charged-Current Interactions of Electrons with Nucleons: This 
requires an additional term on the RHS of equation (||), of the form, 
-5 a ,e8{3, e [Pa,i3(E)r a (E) + (1 - p a ^)T e (E)}, where T a and T e are the appropriate 
differential rates for absorption and emission, from the medium, of a v e . 

D. Neutrino-Neutrino Scattering: At the present stage of the development 
of the formalism the effects of v — v scattering must be put in by hand in a way 
that suits the application. 

In modern supernova calculations in the absence of neutrino mixing, it is usual 
to bypass the explicit Boltzmann equation for the evolution of the density func- 
tions, and to do neutrino transport by following the neutrino distributions nu- 
merically, energy bin by energy bin, and region by region in space. The basic 
input is the differential rate function, T(E, E\) introduced above. Much of the 
review to follow will be concerned with the calculation of these functions in the 
presence of strong interactions among the constituents of the matter. 

However, we note that the equation governing the evolution of the distribution 
function in the presence of mixing cannot be interpreted as assigning probabilities, 
effective immediately, at the outcome of every neutrino interaction. It is clear that 
the numerical simulation of a case with neutrino mixing, involving some phase 
information that may propagate through many scattering interactions, will be 
more complex, and will require going back to the equation for the distribution 
function that is shown above. 



14 Prakash, Lattimer, Sawyer, and Volkas 

3.2 Determination of the Rate Functions 

In applications we need only terms of order G 2 ?. The rates, in media that are in 
thermal equilibrium are given by thermal averages of Heisenberg picture current 
operators. For neutral current scattering we define a current correlator, W^ u , 

W tlu (q,uj)=Z- 1 J d A xe-^e 1 ^ Tr[ e -^-^^(x, t)j v (0, 0)\ , (7) 

where the jn is the neutral current operator, which we have taken to be indepen- 
dent of v flavor .^] The quantities H and iVj contain only the coordinates of the 
species that are in equilibrium. The differential rate for neutrino scattering, at 
angle 6 and with energy loss of u>, or gain of (— u), is 

, f r fl = (^ 2 r l G 2 F (Ei - cv) 2 [l - ME, - u,)]A» v (q, u)W^(q, u) , (8) 
auj a cos 

where q = Pi — P2 is the momentum transferred from the neutrino to the medium, 
Z is the partition function, and 

A"" = (4E.E2)- 1 Tr[^ (1 - 75)7^2 7"(1 " 7s)]- (9) 

For the case of charged current reactions there is a parallel construction, where 
we now take j M to be the hadronic current that is coupled to electron emission, 
with its Hermitean conjugate the current that is coupled to electron absorption: 

W$\q,u>) = -iZ- 1 J d 4 x e -<q-* e i(«+/i)t TV[e-^- s ^^)j M (x,t)it( ,0)] , 

(10) 

where ft is the chemical potential difference jl = /% — fj, p . In the early universe 

application, where the medium is comprised almost entirely of e 1 * 1 , u, u, it is 
5 In applications in which (v e , e) scattering plays a role, it is convenient to include the charged 

current contribution to the [y e , e) interaction in this term, giving a dependence on the v flavor 

indices that we do not show explicitly. 
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quite a good approximation to replace the formal expressions for the functions 

that determine the rates, such as W^ u in equation (||), by their vacuum Born 

approximations, modified only by the statistical factors for the external particles. 

But there are corrections of order 1% in the predicted 4 He abundance that come 

from careful nonequilibrium calculation of the functions related to Wnv i n the 

thermal environment 11 ) . 



In the supernova environment, however, the densities and temperatures from 
the core to the neutrinosphere, at the average radius of the last v scattering, are 
such that the wavelength of a thermal neutrino is greater than the interparti- 
cle spacing, so that neutrino scattering is a collective phenomenon. In this case 
a treatment must begin directly from the dynamical functions defined in equa- 
tions (^) and (|l0|) rather than from the (cross section x density) expression that 
applies to dilute systems or for more energetic neutrinos. In the center of the 
region the densities are on the order of nuclear densities (n Q ~ 0.16 fm -3 ) and the 
matter is quite degenerate, much like the nuclear matter that has been discussed 
with respect to the interiors of heavy nuclei or neutron stars. But nearly all the 
way to the z^-sphere, at a density ~ O.Olno and with T ~5 MeV, strong interac- 
tions among the nucleons will significantly affect the rates of neutrino processes, 
due to the fact that T is less than the per-nucleon interaction energy throughout 
most of the region. 

3.2.1 Long- Wavelength Limits 

The results of the detailed calculations that will be summarized below indicate 
that the main features of the combined effects of Fermi statistics and of the 
nuclear interactions can be described by some simple limits. When the neutrino 
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has energy of the order T <C Mjy there is relatively small energy transferable 

to the nucleons. Then integrating equation (||) over uj to get the rate that is 

differential in angle only, the correlator W^q, u) has a peak around u = and 

the multiplying factors can be evaluated at u = 0. Thus what enters the total 

rate is the the energy integral of the correlator itself, or the t = value of the 

time dependent correlator. As illustration we take the vector current only, with 

a coupling to a single kind of nucleon, and define a structure factor, S(q), by 

/oo 
du W 00 {q,oj) . (11) 
-oo 

Looking first at a free nearly-degenerate Fermi gas, T <C Ep, we find that S(q) 
is less than unity because of the Pauli principle reduction of phase space. There 
is some phase space available for two reasons: 

(a) The Fermi surface is diffuse, since with some states below the Fermi 
energy unoccupied, and some above occupied. No matter how small the mo- 
mentum transfers, this gives some room for transitions, with a rate contribution 
proportional to T. 

(b) There is also a T = limit, but here the neutrino scattering, with momentum 
transfer q, must excite a state from below to above the Fermi level. 

At low temperatures these corrections are additive. But when q ~ T the (a) 
terms, proportional to T, will always predominate over the (b) terms, by the 
ratio of the speed of light to the Fermi velocity (|l2|), (The nucleons are quite 
non-relativistic in most cases.) 

Turning to the case with nucleon- nucleon interactions, again for thermal q, so 
that the terms proportional to T dominate, we can use a powerful classical result: 
in our one species gas the long wavelength limit of the structure factor is given 



Neutrino Propagation in Dense Astrophysical Systems 17 
by 

Lim g _^o S(q) = uTk^ 1 . (12) 

Here kt is the isothermal bulk modulus, kt = n[dP/dn]T = nd 2 F/dn 2 , where P 
is the pressure and F is the Helmholtz free energy density. Note that if the equa- 
tion of state (EOS) were that of a free Boltzmann gas P = nT, equation (^) gives 
the uncorrelated value 5 = 1. If the EOS is that of a degenerate Fermi gas, then 
we get the same answer as in the dominating term proportional to T in the direct 
Pauli blocking calculation reviewed above. Using equation (|i~2| ) we have replaced 
a correlation function calculation by an EOS calculation. This pays dividends for 
the case of the strongly interacting gas. For example, in the nearly degenerate 
case we get the dominating term in the structure function, proportional to T, 
from the zero-temperature EOS, for which nuclear matter calculations provide 
data. We also get a qualitative lesson directly from equation ([12]), namely that 
attractive interactions, which soften the equation of state (reduce kt), increase 
neutrino scattering while repulsive interactions decrease neutrino scattering. 

For the neutrino scattering case, the axial current (Gamow- Teller or GT for 
short) terms coming from equation (^) depend on the nucleon spin density corre- 
lators, rather then the density correlators. These terms provide roughly 3/4 of the 
scattering, in the free particles case. It is straightforward to generalize the above 
small q analysis. The results are that the GT terms in the total opacity have a 
structure similar to equation (|i~2|), but with kt in the denominator replaced by 
K spin = n s d 2 F/dn 2 , where n s is the spin density of the system. That is to say, in 
order to determine the spin density correlators that determine the rate functions 
we must calculate the free energy (or energy, in the case of the nearly degenerate 
system) with a constraint of non-vanishing expectation of nuclear spin, to second 
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order in the spin excess. This would seem to be within the range of variational 

nuclear matter calculations (0), but it has not been carried through for the phe- 

nomenological potentials that are in use. For the case of the nuclear matter in 

nuclei however, to the extent that the term "nuclear matter" is applicable, it is 

the modulus K sp ; n that determines the parameters of the giant GT resonances, 

so that there is indeed the possibility of using values determined fairly directly 

from data rather than those coming from nuclear matter calculations based on a 

specific potential. 

3.2.2 Strong and Electromagnetic Correlations 

With these remarks as orientation, we describe the ring graph calculation of the 
correlators. Although we will quote results from more than one source, and more 
than one exact set of assumptions, the basic inputs are Fermi liquid parameters 
that incorporate a good deal of nuclear phenomenology. These parameters are 
effectively the coupling constants of zero range N-N potentials that describe the 
low energy excitations of the matter. The ring graph sum has two properties that 
are important to us in the light of the above discussion: 

(1) In the q = limit, the sum of the ring graphs for the correlators gives back 
the connection equation (12) when we take the system's free energy density to 



be given by its Hartree value. There is a similar connection for the spin response 
part. 

(2) In contrast to the model of a free nucleon gas, in which the neutrino-nucleon 
collisions are nearly elastic, the ring results show that significant gains and losses 
of neutrino energy may be realized in the collision. Much of this inelasticity 
is found to be in the form of the emission and absorption of spin waves in 
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the medium: these are the continuous matter analogues of the giant GT res- 
onances in nuclei. The giant GT resonance was a principal ingredient in the 
phenomenological determination of the "spin-isospin" Fermi liquid parameter 
that goes into our model. The dispersion relation that describes the nuclear 
GT resonances (with "momentum" ~ 7r/nuclear radius) is consistent with the 



Cerenkov angle found in (14, 15, 1£, 17, 18) for the emission of a spin wave into 



the medium. Whatever residual deficiencies there may be in the method, or in 
our underlying simplification of the nuclear interactions, what follows captures 
some of the essential physics. The ring equations, which we shall give in rather 
schematic form, begin with the definition of a "polarization" IJ^ through the fol- 
lowing replacements on the LHS and RHS of equation (||) W fll/ (q, lo) — > 11^ (q, a;) 
and j^(x, t)j u (0, 0) — ► [j^(x, t),j u (0, 0)]6(t), the retarded commutator being the 
structure of choice for many-body graph summing. Then the rate function W^ u 
is recaptured through, 

WW(?,w) = 2(1 - e-^r^mtrV^,^]. (13) 

The function IJ fJjU (q, u) can be decomposed in scalar functions. For simplicity 
we illustrate the ring approximation (sometimes referred to as the RPA) for the 
case of the vector current part of the neutral current coupling to a single species 
of nucleon (or, equivalently, to the isoscalar part of the coupling to protons and 
neutrons). In this case the time components Uoo(q,oj) = U(q,Lu) obey the ring 
equation, 

= ijgfe - (14) 

where v(q) = J d?x exp[iq • x]F(x) and II^ (q, uj) is the free polarization part. 
If we take only the numerator in equation (14), we recover the effects of Pauli 
blocking alone. 
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This approach can be fleshed out to include spin and isospin, giving rise to 

multichannel equations; then the minimal v(q) in equation (|l4|), now a matrix in 

channel space, can be taken from the Fermi liquid parameters as estimated, e.g., 

in (j]~9|) . Details and further discussion of these methods can be found in Refs. 



The extension of these methods to encompass a relativistic field-theoretic de- 



scription of the matter in which the neutrinos propogate is straightforward (16, 
|l7|,|20|); in this approach the target particle time ordered or causal polarization 
tensor is calculated using 



n Mi , = -i 



I (2^ Tr [ r ( G 2(p)^G 4 (p + g)J,)]. (15) 

The Greens' functions Gi(p) (the index i labels particle species) describe the 
propagation of baryons at finite density and temperature. The current operator 
J a is la f° r the vector current and 7^75 for the axial current. Effects of strong and 
electromagnetic correlations may be calculated by utilizing the RPA polarization 
tensor 

U RPA = u + U RPA DU ) 
where D denotes the interaction matrix (see (|17|,|2"0"D for more details). 

3.2.3 Neutrino Scattering in Heterogeneous Media 

The coherent scattering of neutrinos from heterogeneous media was first dis- 
cussed by Freedman (^) in conjunction with heavy nuclei in the subnuclear 
density regime. For neutrino de Broglie wavelengths X u >> Rn, where Rn is 
nuclear radius, the opacity is enhanced by a factor N 2 /6A rj 100 compared to 
that from the same mass density of nucleons alone (pl[). In the supernova en- 
vironment, coherent scattering dominates all other opacity sources. However, 
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nuclei are relatively close together and the neutrino wavelengths are not neces- 
sarily large. Three important corrections (22) that must be applied which reduce 



the coherent cross section include the nuclear structure factor (when X u < Rn), 
liquid structure effects (when \ v > Rn (no/n) 1 / 3 , the internuclear spacing), and 
polarization of the medium by the electrons (when X u > Rry, where Rry is the 
Debye radius). Obviously, effects that influence the nuclear radius, such as the 
internuclear spacing, the nuclear surface energy, and finite temperatures, are very 
important, and are equation of state sensitive. In addition, around no/3, where 
the nuclear phase might transform into an inside-out, or bubble, phase, structural 
effects could play a role in neutrino cross sections. 

More recently, coherent scattering from droplets of exotic matter (kaon conden- 
sates or quark matter) in the supranuclear density regime has been also considered 
( p3| ) for the PNS environment. 

The Lagrangian that describes the neutral current coupling of neutrinos to the 
droplet (either a nucleus or a droplet of exotic matter) is 

£w = -^j= ^7m(! - 75)^ Jd > ( 17 ) 

where Jj^ is the neutral current carried by the droplet. For non-relativistic 
droplets, J^, = pw{ x ) has only a time like component. Here, pw{%) is 
the excess weak charge density in the droplet. The total weak charge enclosed 
in a droplet of radius Rn is N\y = J^ N d 3 x pw{%) and the form factor is 
F(q) = (1/Nyy) J^ N d 3 x pw{x) smqx/qx. The differential cross section for neu- 
trinos scattering from an isolated droplet is then 

da F 1 ^ 

+««)*«(,). (18) 
where E v is the neutrino energy, 6 is the scattering angle, and q = y/2E w (l— cos 9) 
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is the momentum transfer. Since the droplets are massive, elastic scattering 
dominates. 

The droplet radius Rn and the inter-droplet spacing are determined by the 
balance between surface and Coulomb energies. In the Wigner-Seitz approxima- 



Except for one aspect, we will neglect coherent scattering from more than one 
droplet. If the droplets form a lattice, Bragg scattering will dominate and our 
description would not be valid. But for low density and a liquid phase, inter- 
ference from multiple droplets affects scattering only at long wavelengths. If the 
temperature is not small compared to the melting temperature, the droplet phase 
will be a liquid and interference from scattering off different droplets are small 
for neutrino energies E v > (l/R\y)- However, multiple droplet scattering cannot 
be neglected for E v <^ l/Ryy- The effects of other droplets is to cancel scattering 
in the forward direction, because the interference is destructive except at exactly 
zero degrees, where it produces a change in the index of refraction of the medium. 
These effects are usually incorporated by multiplying the differential cross section 
equation (|l8|) by the static form factor of the medium 



where g(r) is the radial distribution function of the droplets. The droplet cor- 
relations, which determine g(r), are measured in terms of the dimensionless 
Coulomb number V = Z 2 e 2 / ' {SirRy^kT). Due to the long-range character of 
the Coulomb force, the role of screening and the finite droplet size, g(r) cannot 
be computed analytically. We use a simple ansatz for the radial distribution func- 
tion g(r < Rw) = and g[r > Rw) = 1- For this choice, S(q) is independent 
of r. Monte Carlo calculations ( p4[ ) of a simple one component plasma indicate 



tion, the cell radius is Rw = (3/AttNd) 



/ 3 , where the droplet density is Njj. 




(19) 
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that this choice of S(q) is adequate for neutrino energies of interest. 

The simple ansatz for g(r) is equivalent to subtracting, from the weak charge 

density p\y, a uniform density which has the same total weak charge Nw as the 

matter in the Wigner-Seitz cell. Thus, effects due to S(q) may be incorporated 

by replacing the form factor F(q) by 

F(q) - F(q) = F(q) - 3 ~ (qRw) cos qR w 

{qRwr 

The neutrino-droplet differential cross section per unit volume is then 

' ' ,(7 N D ^G 2 F N^(1 + cos 9)F 2 (q). (21) 



VdcosO 16vr 

Note that even for small Njj, the factor Nyy, which in the case of nuclei is 
proportional to N, serves to enhance the droplet scattering. 

This concludes our overview of some of the technical tools that are needed and 
available to address media-related issues. 

3.3 Neutrino Mean Free Paths: Examples 

The differential cross section in equation is required in multi-energy group 
neutrino transport calculations (see Section 5.1). However, more approximate 
neutrino transport schemes (such as those described in Section 5.2.2) utilize the 
total cross section per unit volume (or equivalently the inverse mean free path), 
integrated over the angle 9 and energy transfer uj in equation (JsJ) , as a function of 
the neutrino energy. Examples of such neutrino scattering (that is common to all 
neutrino species) and absorption mean free paths for conditions relevant to the 
deleptonization and cooling epochs of PNSs are considered below. The equations 
of state used in these calculations are described in (f25|). 
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A. Effects of Composition: Under degenerate conditions even modest changes 

to the composition significantly alter the neutrino mean free paths. In Figure |l] 
the top (bottom) panels show the scattering (scattering plus absorption) mean 
free paths in neutrino-free (neutrino-trapped) matter. The left (right) panels 
show results for selected temperatures (for the neutrino energy E u = irT) in mat- 
ter without (with) hyperons. The presence of hyperons significantly decreases 
the mean free paths in both environments because of the additional available 
channels. 

During the deleptonization stage, charged current reactions in fact dominate 
scattering reactions. At zero temperature, reactions like the direct Urea process 
v + n <-» e + p depend sensitively on the proton fraction Y p (|26|). Kinematic 
restrictions require Y p to be larger than 11 — 14% (this is called the direct Urea 
threshold). At early times, when large numbers of trapped neutrinos are present, 
these reactions proceed without hindrance. After several tens of seconds, however, 
Y p , which depends sensitively on the density dependence of the nuclear symmetry 
energy, decreases to its cold, catalyzed value. In field-theoretical models, the 
symmetry energy is largely controled by the p-meson exchange which increases 
strongly with density, and establishes a typical critical density n c = 2 ~ 3uq. 
However, at finite temperature, the equilbrium Y p and average E v values are 
larger than at zero temperature, enabling the charged current reactions to proceed 
even below n c . Figure |2| shows that this is indeed the case even at relatively low 
temperatures (T ~ 3 — 5) MeV for a baryon density = 0.15 fm -3 . Thus, Urea 
processes dominate the opacity until very late times. 
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B. Effects of Strong and Electromagnetic Correlations: The RPA and 

Hartree scattering mean free paths of thermal neutrinos in neutrino-free matter 
are compared in Figure ||. Correlations are generally more important with in- 
creasing density and decreasing temperature. The density dependence is chiefly 
due to variations in the effective baryon mass, which is controlled by the Fermi 
Liquid parameter g' . The RPA correlations, even for a large g' , will typically 
increase the mean free path by a factor 2 — 2.5. 

The RPA and Hartree absorption mean free paths are compared in Figure ^| 
for neutrino-trapped matter. The qualitative trends for neutrino- free matter are 
similar, although the largest enhancements due to correlations are a factor of 
5 in the neutrino-free case instead of about 2 in the neutrino-trapped case. In 
both cases, the results are quite similar to those for the scattering case shown in 
Figure |3|. 

C. Mean Free Paths in Quark Matter: An intriguing possibility is that 
sufficiently dense matter contains deconfined quarks. It is likely that quark mat- 
ter, if it exists, will comprise one component of a mixed phase with hadrons (]27|). 
Within the mixed phase, the thermodynamic and chemical potentials will be 
rather different than in ordinary hadronic matter of the same density. As a 
result, significant effects on neutrino opacities are anticipated. Figure || shows 
scattering and absorption mean free paths, for the cases of matter with trapped 
neutrinos and neutrino- free matter, in which a mixed phase occurs pSj). The 
vertical dashed lines indicate the central densities of 1.4 M stars and maximum 
mass stars, respectively; densities above the right-most vertical line cannot be 
achieved in any neutron star. The thin curves show the mean free paths for the 
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respective matter (hadrons or quarks) in their pure phases, while their values 

within the mixed phase are indicated by thick curves. The chief consequence 
of a mixed phase is that it reverses the trend of pure hadronic matter to have 
steadily decreasing mean free paths with increasing densities. In other words, 
matter within the mixed phase region is substantially more transparent than 
matter without such a transition, whether or not the matter is composed of pure 
hadrons or quarks. Partly, this result is a consequence of the fact that the en- 
tropy is held fixed in these panels: quark matter has a higher specific heat than 
hadronic matter and favors a lower temperature. In the presence of a mixed 
phase, the temperature is smaller than in its absence, and the neutrino cross 
sections, which increase roughly as T 2 , decrease. 

D. Droplet Scattering In Heterogeneous Phases: Models of first order 
phase transitions in dense matter ( p9| ) provide the weak charge and form factors 
of the droplets and permit the evaluation of i/-droplet scattering contributions to 
the opacity of the mixed phase ( p3| ) . For first order kaon condensate and quark- 
hadron phase transitions, the neutrino mean free paths in the mixed phase are 
shown in the left and right panels of Figure [6|, respectively. The transport cross 
sections in diffusive transport are usually differential cross sections weighted by 
the angular factor (1— cos 9). The kaon droplets are characterized by radii Rn ~ 7 
fm and inter-droplet spacings Rw ~ 20 fm, and enclose a net weak vector charge 
N\y ~ 700. The quark droplets are characterized by Rn ~ 5 fm and Rw ~ 11 
fm, and an enclosed weak charge Nw ~ 850. For comparison, the neutrino mean 
free paths in uniform neutron matter at the same density and temperature are 
also shown. A large magnification in the opacity exists for E v ~ ttT. At much 
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lower energies, the inter-droplet correlations tend to screen the weak charge of 
the droplet, and at higher energies the coherence is attenuated by the droplet 
form factor. Similar effects occur in the case of coherent scattering from nuclei 



4 NEUTRINOS IN THE EARLY UNIVERSE 

4-1 Big Bang Nucleosynthesis 

An interesting application of neutrino transport in a dense medium lies in the 
study of neutrino oscillations in that phase of the early universe immediately 
prior to and during Big Bang Nucleosynthesis (BBN) (|30|). Neutrinos play two 
very important roles during BBN. 

First, light neutrinos will be relativistic during the nucleosynthetic epoch, and 
so contribute to driving the expansion of the universe. Under radiation domi- 
nance, the expansion rate or Hubble parameter is given by 



where natural units ft = c = 1 have been employed, T is temperature and Mp is 
the Planck mass. The effective number of massless degrees of freedom is 



where gi counts the internal states of species i. Light neutrinos contribute to the 
sum over fermionic species. If particle of type i is in thermal contact with the 
plasma, then Tj = T. If the expansion rate is smaller than the rate for a given 
process (e.g. a nuclear reaction), then that process will be dynamically important. 
The second BBN role arises from electron neutrinos and antineutrinos trans- 



AD- 




(22) 




(23) 



28 Prakash, Lattimer, Sawyer, and Volkas 

forming protons into neutrons and vice-versa through the reactions 

v e n <-» e~p, V e p <-> e + n . (24) 

When chemical equilibrium obtains, the ratio of neutron to proton number den- 
sities is given by 

^ = e* P (-^p-^), (25) 

where fi Ve is the electron neutrino chemical potential. A crucial event just prior 
to BBN is weak freeze out, when the reaction rates for equation (pj) drop below 
TC. These reactions are then no longer rapid enough to maintain the n/p ratio 
at its equilibrium value, and neutron decay becomes important. The n/p ratio 
determines the amount of primordial 4 He synthesized, because to a good first ap- 
proximation all available neutrons are incorporated into this species. "Standard 
BBN" sees jjL Ue arbitrarily set to zero. However, there is no known fundamen- 
tal reason to require fr Ue = 0. Neutrino oscillations can lead to time dependent 
neutrino chemical potentials, with a direct impact on BBN through the n/p ratio. 

The mean total collision rate {T a ) for neutrino v a (a = e, [i, r) in the epoch of 
the early universe after u/]2 disappearence is 

(r„J ~ Ka G 2 F T 5 , (26) 

where K e ~ 4, k^ iT ~ 2.9 and Gp is the Fermi constant. The average is taken 
over a Fermi-Dirac (FD) distribution 

A-«(P,,,„,T)^/«(f,^) = ^ 7E ^ T , (27) 

where p ~ E is the magnitude of the three-momentum of the neutrino. The 
chemical potential [x Va has been arbitrarily equated to zero in computing the 
righthand side of equation (p6|). If antineutrinos are zero chemical potential FD 



Neutrino Propagation in Dense Astrophysical Systems 29 
distributed, then their mean collision rates are also given by equation (26). Note 

that the total equilibrium number density n Va is given by 

nv a (Vu a ,T)= / N^(p,fx Ua ,T)dp. (28) 
Jo 

If we take the constituents of the plasma at this epoch to be free of exotic states 
such as light sterile neutrinos, light neutral spin-0 bosons and the like, then from 
equations (|22|) and (|23|) the expansion rate of the universe is 

H ^ 5.44— , (29) 
M P v ' 

with = 10.75 under the stated assumptions. 

Neutrino collisions will be very important for all temperatures above the de- 
coupling temperature, T Va , defined through 

K a G 2 F Tl ~ 5.44^J T Va ~ 1 MeV . (30) 

Above about 1 MeV, collisions keep neutrinos in thermal equilibrium with the 
electrons, positrons, photons and the other neutrinos and antineutrinos. Also, 
neutrino oscillations will be collision affected in this regime, necessitating the 
use of the Quantum Kinetic Equation formalism (f|,[5, 31, 32, 33, 34, 35, 36|) . Since 



BBN begins at about 0.4 MeV, the neutrinos are expected to be propagating in 
a collisionless way during the BBN epoch itself. However, the initial conditions 
for BBN (such as the n/p ratio at freeze-out) depend on the collision-affected 
neutrino dynamics of the previous epoch. 

Naively, one would hope to focus on oscillation dynamics involving the three 
known neutrinos v e ,u,T- However, most of the interesting work has focussed on 
active-sterile neutrino oscillations. The reasons for this are, first, that sterile 
neutrinos can have dramatic consequences for BBN, and, second, the Quantum 
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Kinetic Equations for an active-sterile system, equation 



® 



are much simpler 



than for its active-active counterpart, equation (||). 



Active-active oscillations are often neglected in the literature because of the 
common assumption that all neutrino chemical potentials are zero, or at least 
extremely small. If the chemical potentials are in fact nonzero, then an important 
consequence is the inequality of neutrino and antineutrino number densities of a 
given flavor. Define the a-flavor neutrino asymmetry L Ua by 



When thermal equilibrium holds this evaluates to 



where £ Va = fi Ua /T, £ Ua = ^v a /T are dimensionless chemical potentials. This 
equation is exact if £ v = — £ Va , otherwise it is a good approximation for |£| ~ 1. 
It is convenient to scale the neutrino number density with the photon number 
density, and to consider the £'s in place of the /i's, because the scaled quantities do 
not red-shift. While processes such as u a V a «-> e + e~ <-> 77 are occurring quickly 
relative to the expansion rate, the neutrinos and antineutrinos are in chemical 
equilibrium which requires £„ Q + £ u = £ e - + £ e+ = 0, so equation fl32|) is then 



The assumption that v/V chemical potentials are zero, or perhaps of the tiny 
magnitude motivated by the baryon/electron asymmetries, is one of the defining 
features of standard BBN. However, the equal number densities so implied for all 
active neutrino species render neutrino oscillations cosmologically uninteresting, 
simply because equal distributions would just be exchanged through the oscil- 
lations. We will therefore concentrate on cases featuring unequal distribution 




1 



(32) 



exact. 



Neutrino Propagation in Dense Astrophysical Systems 31 
functions. A few well motivated situations of this sort are: 

1. Sterile neutrinos. They would decouple very early, with subsequent reheat- 
ing processes greatly diluting the putative sterile component of the plasma, 
so that by the epoch immediately prior to BBN their number density would 
be effectively zero. Active-sterile oscillations could then repopulate the 
plasma with sterile neutrinos. 

2. Mirror neutrinos. They are similar to sterile neutrinos, except that they 
have weak-scale self-interactions, and their small initial number densities 



are justified differently @|3|,||,gg) 



3. Active neutrinos in the presence of large chemical potential differences be- 
tween the flavors. As we will see, active-sterile oscillations can create rea- 
sonably large chemical potentials during the epoch preceeding BBN (|4l|). 
After this has happened, other oscillation modes, including active-active 



channels, can reprocess the flavor of the lepton asymmetry (42). In addi- 



tion, large chemical potentials can be created during much earlier epochs, 



for example by an Affleck-Dine mechanism 44), which can be similarly 
reprocessed. 

4-2 Active- Sterile Oscillations: Formalism 
4.2.1 Quantum Kinetic Equations 

Equations (||) and (||) describe the evolution of the 1-body reduced density matrix 
p for an active-sterile system. The antineutrino system is described by a similar 
equation for p. 

The diagonal entries of p (p) are appropriately normalized distribution func- 
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tions for v a (V a ) and v s (V s ): 

N Ua (P,T) = ~[P (y,T)+P z (y,T)}N e<! (p,0,T), (33) 
N„.(p t T) = ~[P (y,T)-P z (y,T)]N e i(p,0,T), (34) 

where N eq (p, 0, T) is the often- used reference distribution function. The antineu- 
trino distribution functions are given similarly. Since we have defined Pq z to 
be ratios of distribution functions, there is no term in the QKEs related to the 
expansion of the universe. It is understood, however, that a neutrino which has 
momentum p at temperature T will red-shift to momentum p' at temperature T 1 
such that p/T = p'/T'. In other words, the scaled momentum 

V = f (35) 

is the natural and most convenient variable. Note that y has nothing to do with 
the unit vector y. 

The function V(y,T), describing the quantally coherent part of the evolution, 
is given by 

V(y,T) = f3(y,T)x + \(y,T)z, (36) 

with 

P(y,T) = —mn29 , \(y,T) = -— cos20 o + V a (y,T) , (37) 

where Am 2 and 9q are, respectively, the mass-squared difference and vacuum 
mixing angle for v a — v s oscillations. The mass eigenstate neutrinos v a ^ are 
defined by v a = cos 0QU a + sin Bo^b, v s = — sin d§v a + cos 9oVb, with 6>o defined so 
that cos 2^o > and Am 2 = m 2 — m 2 . The function V a is the effective matter 



potential (45,|4q). Calculated to the appropriate order for our applications, it is 
given by @ 

\m 2 

V a (y,T) = -—[-a(y,T) + b(y,T)}, (38) 
2yT 
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with the dimensionless functions a(y,T) and b(y,T) being 

■fe^ = - «y < f> ->, *,T) = - mg±gpr_ , (39, 
7r z Z\m/ 7r z Am z myy 

where m^y is the V7-boson mass, A e ~ 17, ~ 4.9 and the a-type effective 

neutrino asymmetry is defined by 

L^^L^+L^ + L^+L^+r,. (40) 

Observe that V depends on p through the dependence of a on L Ua , making the 
asymmetry evolution non-linear. The quantity rj is a small term related to the 
cosmological baryon-antibaryon asymmetry. For antineutrinos, the corresponding 
function V is obtained from V by replacing by —L^ a \ 



The Mikheyev-Smirnov-Wolfenstein (MSW) resonance conditions are (45, 4' 



H 

cos2# + a(y,T) -b(y,T) = 0, cos 20 o - a(y, T) - b(y, T) = , (41) 

for neutrinos and antineutrinos, respectively. It is important to appreciate that 
the resonance conditions at a given temperature are met only for neutrinos (an- 
tineutrinos) of a certain momentum y res (y res ). 

The decoherence function is D(y,T) = T Ua (y,T)/2, with T Ua (y,T) being the 
total collision rate for v^s with momentum y. When thermal and chemical 
equilibria hold, 

r«(y, T) = Ka ^^l G 2 F T 5 y + O(L) , (42) 

where the O(L) term takes care of possible asymmetries in the medium. The 
antineutrino decoherence function is obtained by reversing the signs of all the 
asymmetries. 

The function D is in general a dynamical quantity because the collision rates 
depend on the neutrino distribution functions and hence on p. However, it ap- 
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proximates to the external function displayed above provided that thermal equi- 
librium obtains and all lepton numbers are small (but not necessarily tiny). This 
greatly simplifies the numerical solution of the equations. Fermi factors have been 
neglected in the calculation of D, and also the repopulation function R discussed 
below. 

The repopulation function R is approximately given by 

R(y,T) = r Q (y,T)| ^^|^ - l[P (y,T) + P z (y,T)]} , (43) 

when all distribution functions except that for v a are of FD form, and the u a 
distribution is approximately FD. The physical interpretation of this expression 
is that all weak interaction processes involving v a are tending to send its ac- 
tual distribution function towards equilibrium FD form. For antineutrinos, R is 
obtained from R by replacing £„ a by £ v and the P's by P's. 

4.2.2 The Adiabatic Limit 

The Quantum Kinetic Equations in principle provide a complete solution to os- 
cillating neutrino dynamics in this context. Indeed, numerical solutions to these 



equations have been computed for certain situations (|40] , |48| , [4g| , p0|) . (Numer- 
ical work which considers the simplified case where the full energy spectrum 
of neutrinos is replaced by the mean momentum can be found, for instance, 
in Refs. (^,41,^,52.^, 54). For calculations in a collisionless regime, see 
Refs. (H,|^,||).) Before reviewing some of these results, it will be helpful 
to extract some analytical understanding of important features of the dynamics 
[l],^8|, 5^, 5£, 30). It turns out that taking the adiabatic limit is a useful analyt- 



ical strategy. The functions (3, A, D and R all depend on time or temperature. 
The adiabatic limit is useful when the rates of change of these quantities are 
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sufficiently small flSq). 

Consider first the higher temperature regime where collisions maintain the f a 's 

(and F a 's) close to thermal equilibrium, so that we can adopt the instantaneous 

repopulation approximation, whereby 

P (y,T) + P z (y,T) N^(p^ a ,T) 

2 N e i(p,0,T) { ' 

is maintained at all times, though possibly with a time-dependent chemical po- 
tential. It turns out that the further approximation R ~ is reasonable in this 
regime, although its use is not strictly necessary in developing the adiabatic limit 
(|60|). However, the extra algebra involved when R ^ obscures rather than 
clarifies, so we will set R = for pedagogical reasons. Note also that R ~ 
is a valid approximation for T < 1 MeV because weak collision effects are then 
unimportant (D ~ also holds in this case). 



With R ~ 0, the QKEs simplify to 
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(45) 



or, in a self-evident matrix notation, 



(46) 



We have dropped the explicit (y,T) dependence for notational simplicity. 

To solve equation (E6|), we first introduce the instantaneous diagonal basis 
through 

Q=UP, (47) 
where IA is a time-dependent matrix that diagonalizes /C, 



JCd = diag(A;i, fo, ^3) = UKU 1 



(48) 
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with £1,2,3 being eigenvalues. In the instantaneous diagonal basis, equation ( |46| ) 

becomes 

The adiabatic limit of equation (|46|) is defined by the approximation 



U^- ~0. (50) 

Equation ( [49|) is then formally solved to yield 

P(t) =W- 1 (t)e^o /Cd(t ' )<i *'w(0)7'(0). (51) 

The time dependences of the various quantities have been reinstated for obvious 
reasons. 

Under most circumstances, for instance |A| ^> D and \(3\ 3> D individually, the 
eigenvalue spectrum consists of a complex conjugate pair (56,59), 



fci,2 = — d ± iuj , (52) 



together with a negative (real) eigenvalue, 



fc 3 = ~ 3o I , o • (53) 



d 2 + 



The easily derived relations, 



d = D + (54) 
w 2 = A 2 + /3 2 + fc 3 ^ + ^|, (55) 

should also be noted. It is useful to call d the oscillation-affected decoherence 
function and u the matter- and collision- affected oscillation frequency. The very 
important third eigenvalue k% is the relaxation rate for the v a «-> za, process (see 



below). For many applications is small, so that d ~ D and cj ~ \M 2 + /^ 2 - I n 
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order to use equation (|5l|), the diagonalization matrices must also be computed. 

The explicit expressions can be found in Refs. (56,p9). 



The adiabatic limit solution supplies a useful picture for the dynamics (in its 
domain of applicability (|56|) of course). (Some of the following insights were 



obtained earlier using a related approach (|4l|,|48|).) Consider first the case of 



zero decoherence, D = 0. This situation obtains after neutrino decoupling. The 



eigenvalues k{ are then (|56| , 



h = k* = + A 2 , k 3 = , (56) 

where u> = \J [3 2 + A 2 is in this case exactly equal to the usual matter-affected 
oscillation frequency. The adiabatic solution of equation ([5l]) is then just a re- 
statement of the usual adiabatic matter-affected evolution one obtains by solving 
the Schrodinger Equation. For instance, one can show that equation (51) becomes 

P z (t) = cos20 m (*)cos20 m (O) - (57) 

with the initial conditions P x ,y(0) = and -P 2j o(0) = 1 (i.e. no sterile neutrinos 



initially). The function cos 29 m = A/ \J A 2 + (3 2 is the cosine of twice the standard 
matter-affected mixing angle. The oscillatory terms arising from ki 2 have been 
taken to average to zero. This is the standard adiabatic MSW result. For instance 
if the effective matter potential (and hence A) varies from +00 at t = to —00 
at t = 00 through the resonance A = 0, then P z evolves from 1 to —1, signalling 
complete u a — > v s conversion. A simple physical picture then results: as neu- 
trinos with momentum y cross a narrow MSW resonance centered at y = y res , 
adiabaticity guarantees full v a <-> v s conversion. (Similarly antineutrinos, but 
note that y res / y res when ^ 0.) For a narrow resonance, the QKEs can be 
approximated by equations which relate the rates of change of the distribution 
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functions to the speeds at which the resonance momenta move through the dis- 



tributions (42, 31). In practice, such an approximation can often also be validly 
used for temperatures slightly higher than 1 MeV. 

Now consider the opposite case, where collisional effects dominate. The real 
parts of fci t 2 then completely damp the oscillatory behaviour driven by the imag- 
inary parts, and a nonzero k% enters the game. Equation ( |5T|) then allows one to 
solve for the coherence function P y in terms of P z , yielding 

Py(t)^^P Z (t). (58) 

(One can also solve for P x .) Substituting for P y in equation ([45]), one obtains 

BP 

~0f * hPz (59) 

as a self-contained differential equation for P z and hence for the distribution 
functions. (A similar equation follows for antineutrinos.) Because the coherences 
P X) y have been eliminated, one can speak of the classical Boltzmann limit of the 
QKEs. Because k% is negative, and under most interesting circumstances much 
smaller in magnitude than the decoherence rate D, equation (j59|) describes a 
relatively slow collision-dominated evolution of P z towards zero, that is, towards 
N Ua = N Ug . It can be viewed as a relaxation effect. 

4.2.3 Neutrino Asymmetry Evolution 

Neutrino asymmetries and chemical potentials are important quantities. As pre- 
viously explained, active-active oscillations are expected to have strong effects 
only if e, [x and r asymmetries are unequal, and an e-like asymmetry directly 
affects BBN through the proton-neutron interconversion reactions. Furthermore, 
the effective matter potential of equation d38| ) has a term proportional to a linear 
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combination of asymmetries. Large asymmetries therefore imply small matter- 
affected mixing angles 9 m , and hence suppressed oscillations. This is of particular 
concern for sterile neutrino models: sufficiently large asymmetries will suppress 
active to sterile oscillations with important implications for the consistency of 
these models with BBN. 

Using the QKEs together with ce+s lepton number conservation, one can easily 
show that 

= 2^/ m-PyW q (p^T)dp (60) 

is the evolution equation for L Va under the action of v a <-> v s and V a <-> V s 
oscillations. We will study this simple two-flavor toy model before considering 
more realistic scenarios. 

The QKEs can be numerically integrated to yield asymmetry evolution curves. 
(For numerical reasons, it is useful to include equation ( p0| ) in the coupled system, 
even though it is redundant). Let us suppose that all neutrino asymmetries start 
out being small, perhaps of the order of the observed baryon asymmetry ~ 10 -10 . 
There is no guarantee that this is realistic, but up to now most of the interesting 
work in the literature has focussed on this case. When Am 2 < and the vacuum 
mixing angle is small, there is a large range of parameters for which explosive 
asymmetry growth begins at a certain critical temperature T c , as displayed in 
Figure [7|. This is a runaway positive feedback phenomenon that occurs when an 
— o approximate fixed point changes from being stable to unstable at T = T c 
(41,4^). After a short spurt of quasi-exponential growth, the evolution settles into 
a less dramatic but still significant approximate power law phase, before reaching 
a steady state plateau at a value in the range 0.2 — 0.35 (depending on the 
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oscillation parameter choice) (^2|) .^] The baryon asymmetry in the plasma, which 

couples in through the 77-term in L^ a \ seeds the neutrino asymmetry growth. The 

existence of the asymmetry growth effect must be carefully taken into account 

when calculating sterile neutrino production for BBN purposes. 

The main features of the asymmetry growth curves have been understood, and 

are discussed in some detail in Refs. (41, 42, 48|, 63). In a nutshell, the explo- 



sive growth phase is governed by the collision dominated dynamics leading to 
equations ( |58|) and (|59|), while the approximate power law regime has adiabatic 
MSW transitions as the driving force. It has been shown that asymmetry growth 
resembling that shown in Figure ^ will occur for the parameter space region 

Am 2 < with |Am 2 | ~ 10" 4 eV 2 , 
10- 10 £ sin 2 20 o £ few x lO" 5 (^) 1/2 . (61) 

For I Am 2 1 < 10 -4 eV 2 , asymmetry can be generated, but it tends to be oscillatory 
(|3l | |5l l |55l , |57l). Values of sin 2 20$ smaller than 10 -10 are too small to generate an 
asymmetry, while values exceeding the upper bound above lead to copious sterile 
neutrino production which is also detrimental to asymmetry growth. One can 
also show that the critical temperature T c is roughly given by 

T c ^(16^20MeV)( |Am2 e l v C r 2g ° ) 1/6 - (62) 

4-3 Active-Sterile Oscillations: Applications 
4.3.1 Sterile Neutrino Production 

For the moment, let us continue to study the toy universe where only v a <-> v s 

and V a <-> V s oscillations occur (for a particular a). An interesting question is: 
6 Reference (|62j) found much smaller values for the final asymmetries. A critique of this work 



can be found in Ref. (p3h . See also Refs. (KMlpj) 
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for a given choice of Am? and 60, how does sterile neutrino production affect 

BBN? For a = /i or r the effect is restricted to potentially increasing the energy 

density in relativistic species and hence the expansion rate of the universe. (Note 

that for Am 2 ~ 10 -7 eV 2 significant Un T <-> v s oscillations can occur only after 

neutrino decoupling, so the overall energy density is unaffected.) For a = e, 

the possible effects are more complicated. Because of the reactions in equation 



(24), the possible existence of an e-like asymmetry must be carefully considered, 
and one has to take into account the effect of distortion and depletion of the v e 
spectrum below the decoupling temperature (36. 55, p7). 

The detailed results will clearly depend on whether or not the plasma has 
large neutrino asymmetries. As we have just seen, an a-like asymmetry will be 
generated during the T ~ 10's of MeV epoch by oscillations if the parameters 
obey equation (|H|). Other mechanisms, such as Affleck-Dine, can create large 
asymmetries at much higher temperatures. Unfortunately, prospects for directly 
detecting the relic neutrino sea remain remote, so there is much room for theo- 
retical speculation. 

The important pioneering works on sterile neutrino production were performed 



before the asymmetry generation effect was understood (|6g,|63,|6^,|6^,|7y,|7^). By 
neglecting the baryon asymmetry and setting all neutrino asymmetries to zero, 
one obtains a simplified system driven by identical neutrino and antineutrino 
oscillation dynamics which are decoupled from each other. For Am 2 > 0, this 
simplification is consistent. However, for Am 2 < 0, the L^ a ' = approximate 
fixed point becomes unstable and the zero asymmetry constraint is not justified. 

A useful summary of the pioneering results is contained in Figures 2 and 4 from 
Ref. (p9|). These plots show contours in the Am 2 — sin 2 29$ plane corresponding to 
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different values of g ma x , a parameter which quantifies sterile neutrino contribution 

to the expansion rate, equation (|22"|), through an effective increase in the number 
of fermionic degrees of freedom contributing to g* . (We will call this parameter 
AN e ff from now on. Note that for a = e its meaning is more subtle (|69D.) The 
precise constraint one should use is at this stage unclear, because of uncertainties 
in the primordial light element abundance observations. Nevertheless, these plots 
provide interesting information for neutrino model builders. For instance, a Un — > 
v s solution to the atmospheric neutrino problem with Am 2 ~ 1CP 2 ' 5 eV 2 and 
sin 2 28q ~ 1 would appear to go hand-in-hand with a fully thermalized v s prior to 
BBN. Such a situation would be cosmologically disfavored because of concomitant 
4 He overproduction. 

As already noted, for Am 2 < one has to take into account the asymme- 
try generation that occurs simultaneously with sterile neutrino production. The 
qualitative effect is clear: because asymmetries suppress oscillations, the rigor- 
ous Am 2 < bounds should be weaker than those displayed in Figure 4 of Ref. 
(|69|), for example. However, detailed calculations (|48|,[7^) reveal that the asym- 
metry amplification process does not modify the "traditional" bounds by very 
much, except when the BBN constraint on AN e ff is taken to be rather stringent, 
say AN e ff < 0.1. The reason for this is that sterile neutrino production tends 



to delay the onset of asymmetry growth (see Figure 12 of Ref. (72) for a clear 
illustration) . 

All of the above related to the artificial two-flavor u a /v s case. But the most 
dramatic consequence of asymmetry creation occurs in more realistic multiflavor 
situations. This is simply because an asymmetry created by one active-sterile 
oscillation mode can suppress a different active-sterile mode (El, 48, 4£, 72). The 
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most studied case was motivated by the atmospheric neutrino problem. As noted 

above, the Vn — * v s solution to the atmospheric neutrino problem appears at first 

sight to be cosmologically disfavored. However, tiny mixing between the v s and 

a more massive v T can completely change the dynamics, because the nonzero L Vt 

created by the small-angle, negative Am 2 v T — ► v s mode can suppress the z/« — > v s 

oscillations that would be very strong otherwise. Explicit calculations show that 

if the mass eigenvalues associated with the v^jvg subsystem are much less than 

an eV, then a v T mass at the few eV level or higher is required (41. 18. 4S, 72). 

(Though somewhat large values were obtained in Ref. (|73|) , the later independent 

results of Ref. ( [72] ) agree with the original results of Refs. (41,|4S|).) 

4.3.2 e-Like Asymmetry and 4 He Abundance. 

Suppose a small-angle active-sterile mode with a relatively large but negative 
Am 2 creates a large asymmetry of a certain flavor fairly early. Later on, it 
is certainly possible for other oscillation modes to reprocess the flavor of the 
asymmetry. For instance, small-angle v T — > v s oscillations might create a large 
L Vt ~ 0.3 which gets reprocessed into a reasonably large (say ~ 0.01) L Ue by 
v T *-* v e oscillations. Recall that the usually innocuous active- active modes can 
become important after flavor asymmetries get created. 



Several scenarios of this type have been investigated in the literature (40, 42] , 61 
One has to study the peculiarities of each neutrino model on a case-by-case 
basis. The more flavors one has, the more complicated the analysis becomes. The 
most ambitious attempt to date relates to the mirror matter model with three 
active and three mirror neutrino flavors (f40|). 



By way of example, Ref. (42) considered a model with the three active neu- 
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trinos and one light sterile flavor with the mass hierarchy m Vr 3> m v , m Ve , m Vs . 

A r-flavor asymmetry is first created by the v T <-> v s mode, and subsequently 
reprocessed into e- and //-like asymmetries by v T <-> v e and v T <-> oscilla- 
tions. The e-like asymmetry is of most interest because of its role in the neutron- 
proton interconversion reactions. Final values for L Ve around 0.02(T I/e /T 7 ) 3 for 
AM 2 > 10 eV 2 were found, where AM 2 is the squared-mass splitting between 
v T and the lighter flavors. The effect of such an asymmetry on 4 He abundance 
is roughly equivalent to reducing the expansion rate by, effectively, half of a neu- 
trino: AN e ff ~ —0.5. Note that the sign of the asymmetry is crucial here: 
positive L v Js reduce the 4 He yield (equivalent to reducing the expansion rate), 
while negative values increase the yield (equivalent to increasing the expansion 
rate). The sign of L Ue is controlled by the sign of the L Ut originally created, which 
in turn is controlled by the unknown initial (high T) values of the asymmetries. 
It turns out that the results quoted above assume initial conditions that lead to 
a positive L Ut . 

Clearly, the relevance or otherwise of these studies for understanding nature 
hinges on the as yet unclear existential status of light sterile neutrinos. Future 
terrestrial experiments such as SNO, MiniBOONE and the long baseline projects 
will help to clarify the situation. Better information about the mass and mixing 
angle spectrum for neutrinos will be needed, as well as better information about 
primordial light element abundances. Precision cosmic microwave background 
measurements also have a role to play by helping to pin down the baryon density 
(a crucial input into BBN calculations), as well as the expansion rate of the 



universe at the time of photon decoupling (75,76, 77). 
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4-4 Remarks About Active-Active Oscillations 

Collision-affected active-active oscillations have received very little attention be- 
cause the relevant Quantum Kinetic Equations are rather complicated, as can be 
seen by comparing equation (|j) with equation (|5|). The studies that have been 
done typically focussed on collision- unaffected adiabatic transitions, for which a 
simpler treatment is possible. The works alluded to in the previous subsubsec- 
tion in fact used such an approach. The collision-affected case remains largely 
virgin territory: the required QKEs have been written down (|],||), but remain 
unsolved, except for a certain interesting special case (|||7||). The special situa- 
tion is that of propagation through a "flavor-blind" medium. Examples include 
an active-active oscillating system traversing a dense nucleonic medium, and a 
fu/fr system propagating through an electron-positron plasma. In both cases, 
the microscopic collisions do not distinguish between the flavors of the oscillating 
system. While the practical relevance of these systems is not entirely clear, a 
very interesting phenomenon has been discovered: synchronization of modes in 
the rapid collision limit. In vacuum, and in a refractive medium, neutrino oscilla- 
tion lengths generally depend on energy. It has been found that rapid flavor-blind 
collisions force all modes to have the same oscillation length, which has been an- 
alytically proven to be a kind of thermal average of the oscillation lengths the 
neutrinos would have in the absence of the collisions. This leads at a sharpening of 
MSW transitions, and it reflects an interesting persistence-of-quantal-coherence 
effect which may find application outside of the neutrino domain (|9|). 
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5 NEUTRINOS IN SUPERNOVAE AND PROTO-NEUTRON 

STARS 

The current gravitational core collapse supernova paradigm is based on the orig- 
inal suggestion by Colgate and White (|79|) that the explosion derives from neu- 
trino energy deposition. However, since that time, much has been learned about 
neutrino-matter interactions that have modified the original model. It has also 
been realized that the implementation of accurate neutrino transport is critical, 
since the deposited energy is a small fraction (< 1%) of the available gravitational 
binding energy of the remnant (~ 3GM 2 /5R « 3 • 10 53 (M/1.4 M Q ) 2 ergs). 

The initial phases of the supernova begin with the destabilization and collapse 
of the core of a massive star (see Ref. ( |80[ ) for the important physics which occurs 
during this period). The collapsing core, which is initially composed of iron peak 
elements with a net electron content Y e = n e /riB — 0.41 — 0.43, divides into 
two portions: an inner, homologous (with infall velocity roughly proportional to 
the radius) core and an outer region that collapses supersonically. The infalling 
matter maintains a nearly constant entropy per baryon s ~ 1; the entropy gain 
from out-of-equilibrium weak interactions is balanced by the energy loss from 
escaping neutrinos. During the collapse, electrons and protons are converted into 
neutrons and neutrinos as the matter attempts to maintain beta equilibrium. 
When the central density reaches about 10 12 g cm -3 , neutrinos are unable to 
escape on dynamical timescales and are essentially frozen, or 'trapped', in the 
matter. The lepton number thereafter remains fixed, at a value of Yl = Y e +Y Ue ~ 
0.4. 

The collapse continues until the central density exceeds no when the increased 
pressure from strong interactions reverses it. A shock is formed at the outer edge 
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of the inner core (where sound waves from the center accumulate) and begins 

to propagate through the infalling outer region of the core. But this shock, 
called the "bounce" shock, stalls at a distance from the center of about 150 km, 
due to the energy expended in nuclear dissociation and neutrino losses. The 
stalled shock becomes an accretion shock which separates supersonically infalling 
matter from hot matter slowly settling onto the inner core, a PNS. After 10-20 
ms, the overall structure evolves quasi- hydrostatically. The nucleons produced 
by the dissociation of heavy nuclei are heated by neutrino absorption from the 
hot, newly formed PNS. Because neutrino emission varies as T 6 , there is a point, 
known as the "gain radius" , at which the heating exceeds the cooling. 

Current calculations differ to some extent as to the outcome of this scenario, 
due to varying input physics and level of approximation to the neutrino trans- 
port problem. In addition, the negative entropy gradient that naturally exists 
is unstable against convection. This convection, which is neutrino-driven, seems 
to eventually assist in reviving the shock in some models. However, convection 
can only be realistically modeled in three dimensions, a task that is only be- 
ginning to be addressed. Additional fluid instabilities may also arise near the 
neutrinosphere, at which the neutrino optical depth to infinity is of order unity 
and where the neutrinos can begin to freely escape the star. The shock posi- 
tion during this hydrostatic epoch is determined by a delicate balance between 
thermal pressure caused by neutrino heating and the ram pressure of infalling 
matter (|8l|). A successful supernova results if this balance becomes unstable, 
which could occur if the accretion shock can be maintained at a sufficient dis- 
tance for a long enough time. The ram pressure decreases as material from less 
dense regions of the outer core is encountered, and could be eventually overcome 
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by the more steady neutrino radiation from the core. Nevertheless, explosions in 

model simulations, even when they occur, appear to be marginal at best. Aside 

from uncertainties stemming from simulating 3-D, general relativisitic, neutrino 

transport, the initial structure of the pre-collapse core, the equation of state, and 

neutrino opacities and emissivities, all play roles in the outcome (js^), 

5. 1 Role of Neutrinos in Gravitational Collapse Supernovae 

Neutrino transport in the supernova environment is described by a Boltzmann 
tranport equation, derivable from the kinetic equation (||) by retaining only the 
diagonal elements in the density matrix p(r,r'). Even with this simplification, 
it is a nonlinear integro-partial differential equation that describes the time rate 
of change of the neutrino distribution function /. Advances made to date in the 
numerical solution of this equation in the supernova context may be found in 
Refs. (|83|). Historically, multigroup methods (in which the equation is discretized 
in energy groups) have involved moment expansions. When the temporal deriva- 
tive of the first moment of / is set to zero, a diffusion equation is obtained, 
but this cannot adequately handle the free-streaming regime at low densities. 
Flux limiting schemes have been used to bridge the diffusive and free-streaming 
regimes, but these are somewhat arbitrary and accurate calibration depends upon 
neutrino opacities and dynamics. In addition, there is a problem with the cou- 
pling of different neutrino-energy groups, especially because of v — e scatterings, 
which involve large energy transfers. An additional complication in supernovae 
is that the approach to thermal and chemical equilibrium, and the conversion of 
diffusive flow to free streaming, occur simultaneously in space and time. Even 
with modern parallel supercomputers, it is necessary to integrate the Boltzmann 
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equation over solid angles to reduce the dimensionality of the problem, with a 

corresponding loss of information about the neutrino angular distribution func- 
tion. This could be important in regimes in which neutrino-driven convection, a 
3-D phenomenon, is occuring. 

Crucial weak interaction processes in the supernova environment include 

p + e-^n + v e , (A,Z) + e~ (A,Z- 1) + i/ e , (63) 

u+(A,Z)^v+ (A, Z) , (64) 
v + e~^u + e~, v+(A,Z)-+v + (A,Z)*, (65) 

e+ + e~ -H/ + P, (A,Z)* -► (A,Z) +v + v, (n,p) -► (n,p) + v + v . (66) 



Reactions (63) begin the process of neutronization and decrease of Yl, whose 
value after trapping determines the masses of the homologous core and initial 
PNS, and thus the available energy for the shock and subsequent neutrino emis- 
sions. The equation of state also influences these quantities, most importantly 
through the nuclear symmetry energy. 

In the subnuclear density regime, the coherent scattering reaction fls3|) from 
nuclei in a lattice is the most important opacity source (see Section 3.2). The 
reactions ( |65|) are important in changing the neutrino energy, and in achieving 
thermodynamic equilibrium. As referred to in Section 5.1, the large energy trans- 
fers of these processes is a hurdle for numerical calculations. The reactions ( |66| ) 
are also important in achieving thermodynamic equilibrium. The bremsstrahlung 
(n + n — ► n + n + u + v) and modified Urea (n+p — ► n + n + e + + v + v) processes 
involving nucleons dominate in many circumstances. For example, the produc- 
tion and thermalization of \x and r neutrinos, which receives contributions from 



all the reactions (66), is dominated by nucleon bremsstrahlung for n > 0.005 and 
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T < 15 MeV The modified Urea process dominates the cooling of PNSs if 

direct Urea processes involving nucleons, hyperons or other strange particles do 

not occur. 

5.2 Neutrinos From Proto-Neutron Stars 

A PNS is born in the aftermath of the gravitational collapse of the core of a 
massive star accompanying a successful supernova explosion. During the first 
tens of seconds of evolution, nearly all (~ 99%) of the remnant's binding energy 
is radiated away in neutrinos of all flavors (85|, 86, 87, 38, 3£, 9C). The neutrino 



luminosities and the emission timescale are controlled by several factors, such as 
the total mass of the PNS and the opacity at supranuclear density, which depends 
on the composition and EOS. One of the chief objectives in modeling PNSs is 
to infer their internal compositions from neutrino signals detected from future 
supernovae like SuperK, SNO and others under consideration, including UNO 



(91). 



5.2.1 General Description of the Birth of Proto-Neutron Stars 

The evolution of a PNS proceeds through several distinct stages (^,0) and with 



various outcomes (25), as shown schematically in Figure ||. Immediately following 
core bounce and the passage of a shock through the outer PNS's mantle, the star 
contains an unshocked, low entropy core of mass ~ 0.7 M Q in which neutrinos 
are trapped (stage 1 in the figure). The core is surrounded by a low density, 
high entropy (5 < s < 10) mantle that is both accreting matter from the outer 
iron core falling through the shock and also rapidly losing energy due to electron 
captures and thermal neutrino emission. The mantle extends up to the shock, 
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which is temporarily stalled about 200 km from the center prior to an eventual 

explosion. 

After a few seconds (stage 2), accretion becomes less important if the supernova 
is successful and the shock has ejected the stellar envelope. Extensive neutrino 
losses and deleptonization will have led to a loss of lepton pressure and the col- 
lapse of the mantle. If enough accretion has occurred, however, the star's mass 
could increase beyond the maximum mass capable of being supported by the hot, 
lepton-rich matter. If this occurs, the remnant collapses to form a black hole and 
its neutrino emission is believed to quickly cease fl92[) . 

Neutrino diffusion deleptonizes the core on time scales of 10-15 s (stage 3). 
Diffusion time scales are proportional to R 2 {c\ v )~ l , where R is the star's radius 
and \ v is the effective neutrino mean free path. This generic relation illustrates 
how both the EOS and the composition influence evolutionary time scales. The 
diffusion of high-energy (200-300 MeV) us from the core to the surface where they 
escape as low-energy (10-20 MeV) vs generates heat (a process akin to joule 
heating). The core's entropy approximately doubles, producing temperatures 
in the range of 30-60 MeV during this time, even as neutrinos continue to be 
prodiguously emitted from the star's effective surface, or u— sphere. 

Strange matter, in the form of hyperons, a Bose condensate, or quark matter, 
suppressed when neutrinos are trapped, could appear at the end of the delep- 
tonization. Its appearance would lead to a decrease in the maximum mass that 
matter is capable of supporting, implying metastability of the neutron star and 
another chance for black hole formation (|25|). This would occur if the PNS's 
mass, which must be less than the maximum mass of hot, lepton-rich matter 
(or else a black hole would already have formed), is greater than the maximum 
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mass of hot, lepton-poor matter. However, if strangeness does not appear, the 

maximum mass instead increases during deleptonization and the appearance of 

a black hole would be unlikely unless accretion in this stage remains significant. 

The PNS is now lepton-poor, but it is still hot. While the star has zero net 
neutrino number, thermally produced neutrino pairs of all flavors dominate the 
emission. The average neutrino energy slowly decreases, and the neutrino mean 
free path increases. After approximately 50 seconds (stage 4), A ~ R, and the 
star finally becomes transparent to neutrinos. Since the threshold density for the 
appearance of strange matter decreases with decreasing temperature, a delayed 
collapse to a black hole is still possible during this epoch. 

Following the onset of neutrino transparency, the core continues to cool by 
neutrino emission, but the star's crust remains warm and cools less quickly. The 
crust is an insulating blanket which prevents the star from coming to complete 
thermal equilibrium and keeps the surface relatively warm (T ~ 3 x 10 6 K) for 
up to 100 years (stage 5). The temperature of the surface after the interior of the 
star becomes isothermal (stage 6) is determined by the rate of neutrino emission 
in the star's core and the composition of the surface. 

5.2.2 The Proto-Neutron Star Evolution Equations 

The equations that govern the transport of energy and lepton number in a PNS 
are obtained from the Boltzmann equation for massless particles (|85, p3, 94, 85). 



We will focus on the non-magnetic, spherically symmetric situation, and note that 
fluid velocities are small enough so that hydrostatic equilibrium is nearly fulfilled. 
Under these conditions, the neutrino transport equations in a stationary metric 

ds 2 = -e 2 ^dt 2 + e 2A dr 2 + r 2 d9 2 + r 2 sin 2 9 d<S> 2 (67) 
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are: 

d{N u /n B ) d(eH<irr 2 F u ) ^Sn 
+ ~ = e. 



dt da ub 

d(Jv/n B ) p d(l/n B ) , ^ dje^Trr 2 ^) _ ^S E 

at at aa 

where n B is the baryon number density and a is the enclosed baryon number 
inside a sphere of radius r. The quantities N u , F u , and Sjy are the number 
density, number flux and number source term, respectively, while J„, ff^, P u , 
and S 1 g are the neutrino energy density, energy flux, pressure, and the energy 
source term, respectively. 

In the absence of accretion, the enclosed baryon number a is a convenient 
Lagrangian variable. The equations to be solved split naturally into a trans- 
port part, which has a strong time dependence, and a structure part, in which 
evolution is much slower. Explicitly, the structure equations are 

— = 1 dm _ p 

da A7rr 2 n B e A ' da n B e^ 
d6 e A / \ dP e A / \ 

8! = 4^(" i + 4 " p ) • *=- ( ' , + p) 4^('" + 4 ' rP )< 71 > 

The quantities m (enclosed gravitational mass), p (mass-energy density), and P 
(pressure) include contributions from the leptons. To obtain the equations em- 
ployed in the transport, equation ( |68|) may be combined with the corresponding 
equation for the electron fraction 

^ = -e^ (72) 
dt n B y J 

to obtain 

«fi + -0. (73) 

dt da 

Similarly, equation (|6^) may be combined with the matter energy equation 

dU +p d(l/n B ) = _ e ,S £ ^ 
dt dt n B 
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where U is the specific internal energy. The first law of thermodynamics yields 



dt 



dt 



da 



(75) 



At high density and for T> 1 MeV, the source terms in the Boltzmann equa- 
tion are sufficiently strong to ensure that neutrinos are in thermal and chemical 
equilibrium with matter. Thus, the neutrino distribution function in these regions 
is both nearly Fermi-Dirac and isotropic. We can approximate the distribution 
function as an expansion in terms of Legendre polynomials to O(p), which is 
known as the diffusion approximation. Explicitly, 



f(u, fi) = /oM + M/lM , fo = [1 + e(^)]- 1 



(76) 



where fo is the Fermi-Dirac distribution function at equilibrium (T = T ma t, 
fj, v = \x e ^\ with u and [i v being the neutrino energy and chemical potential, 
respectively. In the diffusion approximation, /i(/o) becomes d88|) 



h 



-D(oj) 



_-A 9 /0 



dr dr dco 



(77) 



The explicit form of the diffusion coefficient D is given by 



D(u) 



(78) 



The quantity j = j a + j s , where j a is the emissivity and j s is the scattering 
contribution to the source term. The absorptivity is denoted by A a and k\ is the 
scattering contribution to the transport opacity. Substituting 



dr V dr T dr J dui 



(79) 



where r\ = fJ, v /T is the neutrino degeneracy parameter, in equation (|77|), yields 



fi = -D(u)e 



-A 



dr Te't' dr 



dfo 
duj 



(80) 
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Thus, the energy-integrated lepton and energy fluxes are 



F, 



V 



6tt 2 3 dr 

e -A e -^ T 3 r d (T e <t>) 



+ (Te^)D 2 



dr] 
Or 



dr/ 
Or 



(81) 




The coefficients D 2 , D%, and D4 are defined by 



'» = / 

Jo 



'OO 



dx x n D(u)fo(u;)(l - f (u)) , 



(82) 



where x = oj/T. These diffusion coefficients depend only on the microphysics of 
the neutrino-matter interactions. The fluxes appearing in the above equations 
are for one particle species. To include all six neutrino types, we redefine the 



D 2 = D u 2 e + D v 2 £ , D 3 = Dl e - Dl £ , D 4 = + + 4D^ . (83) 



5.2.3 Neutrino Luminosity from Proto-Neutron Stars 

A fair representation of the signal in a terrestrial detector can be found from the 
time dependence of the total neutrino luminosity and average neutrino energy 
together with an assumption of a Fermi-Dirac spectrum with zero chemical po- 
tential. The total neutrino luminosity is globally the time rate of change of the 
star's gravitational mass, and due to energy conservation, is also 



at the edge of the star. However, since the spectrum is not precisely Fermi- 
Dirac at the neutrinosphere, a diffusion scheme only approximates the average 
energy. The average energy can be approximated as < E v >ss 3T U , where T u is 
a mass-averaged temperature in the outermost zone. 




L u = e 2(t> ^r 2 H l 



V 



(84) 
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Neutrino signals from PNSs depends on many stellar properties, including the 

mass; initial entropy, lepton fraction and density profiles; and neutrino opacities. 



In Figures ^-13, the dependence of neutrino emission on PNS characteristics are 
shown from the detailed study of Pons et al. (|88|,|39|,|90|). The generic results (see 



Figure ||) are that both L v and < E v > increase with increasing mass (85, 3^ 



< E v > for all flavors increases during the first 2-5 seconds of evolution, and then 
decreases nearly linearly with time. For times larger than about 10 seconds, and 
prior to the occurrence of neutrino transparency, the L u decays exponentially 
with a time constant that is sensitive to the high-density properties of matter. 
Significant variations in neutrino emission occur beyond 10 seconds: L u is larger 
during this time for stars with smaller radii and with the inclusion of hyperons in 
the matter. Finally, significant regions of the stars appear to become convectively 



unstable during the evolution, as several works have found (95). 

The main effect of the larger mean free paths produced by RPA corrections 
(|T3,[l5,pl is that the inner core deleptonizes more quickly (see Figure [lQj) . In 



turn, the maxima in central temperature and entropy are reached on shorter 
timescales. In addition, the faster increase in thermal pressure in the core slows 
the compression associated with the deleptonization stage, although after 10 s 
the net compressions of all models converge. The relatively large, early, changes 
in the central thermodynamic variables do not, however, translate into similarly 
large effects on observables such as L v and < E v >, relative to the baseline 
simulation. It is especially important that at and below nuclear density, the 
corrections due to correlations are relatively small. Since information from the 
inner core is transmitted only by the neutrinos, the time scale to propagate any 
high density effect to the neutrinosphere is the neutrino diffusion time scale. Since 



Neutrino Propagation in Dense Astrophysical Systems 57 
the neutrinosphere is at a density approximately O.Olno, and large correlation 

corrections occur only above no/3 where nuclei disappear, correlation corrections 

have an effect at the neutrinosphere only after 1.5 s. However, the corrections 

are still very important during the longer-term cooling stage (see Figure |ll]) , and 

result in a more rapid onset of neutrino transparency compared to the Hartree 

results. 

5.2.4 Neutrino Signals in Terrestrial Detectors 

A comparison of the signals observable with different detectors is shown in Fig- 



ure 



12| , which displays L function of baryon mass Mb for stars containing 

quarks in their cores. In the absence of accretion, Mb remains constant during 
the evolution, while the gravitational mass Mq decreases. The two upper shaded 
bands correspond to estimated SN 1987A (50 kpc distance) detection limits with 
KII and 1MB, and the lower bands correspond to estimated detection limits in 
SNO, SuperK, and UNO, for a Galactic supernova (8.5 kpc distance). The detec- 
tion limits have been set to a count rate dN/dt = 0.2 Hz (|89|). It is possible that 
this limit is too conservative and could be lowered with identifiable backgrounds 
and knowledge of the direction of the signal. The width of the bands represents 
the uncertainty in < Ep e > due to the diffusion approximation (88,^9, 90). It 
appears possible to distinguish between stable and metastable stars, since the lu- 
minosities when metastability is reached are always above conservative detection 
limits. 
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5.2.5 Metastable Proto-Neutron Stars 

Proto-neutron stars in which strangeness appears following deleptonization can 
be metastable if their masses are large enough. One interesting diagnostic that 
could shed light on the internal composition of neutron stars would be the abrupt 
cessation of the neutrino signal. This would be in contrast to a normal star of 
similar mass for which the signal continues to fall until it is obscured by the 
background. In Figure the lifetimes for stars containing hyperons (npH), 
kaons (npK) and quarks (npQ) are compared (^). In all cases, the larger the 
mass, the shorter the lifetime. For the kaon and quark PNSs, however, the 
collapse is delayed until the final stage of the Kelvin-Helmholtz epoch, while this 
is not necessarily the case for hyperon-rich stars. In addition, there is a much 
stronger mass dependence of the lifetimes for the hyperon case. 

Clearly, the observation of a single case of metastability, and the determination 
of the metastability time alone, will not necessarily permit one to distinguish 
among the various possibilities. Only if the metastability time is less than 10-15 
s, could one decide on this basis that the star's composition was that of npH 
matter. However, as in the case of SN 1987A, independent estimates of Mb 
might be available (|96|). In addition, the observation of two or more metastable 
neutron stars might permit one to differentiate among these models. 

6 OUTLOOK 

On the early universe front, an important issue is simply whether or not light 
sterile neutrinos exist. If they exist, then the details of light element synthesis 
could be interestingly different from that of standard Big Bang Nucleosynthesis, 
because of the neutrino asymmetry amplification phenomenon. One must look 
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forward to crucial experimental results from SNO, MiniBOONe and the long 

baseline facilities. Naturally, further observational study of the actual light ele- 
ment abundances is a central concern. Precision cosmic microwave background 
anisotropy measurements will also be important as an independent probe of the 
baryon to photon ratio, and of the expansion rate of the early universe. An in- 
teresting theoretical issue currently under examination is that of inhomogenous 



neutrino asymmetry creation (97,98). 

The outlook from the supernova perspective hinges on technical advances in 
handling multigroup, general relativistic, Boltzmann neutrino transport. The 
advent of next-generation neutrino detectors such as Super-Kamiokande and the 
Sudbury Neutrino Observatory promises thousands of neutrino events in the next 
Galactic supernova. These will provide crucial diagnostics for the supernova 
mechanism, important limits on the released binding energy and the remnant 
mass, and critical clues concerning the composition of high density matter. Re- 
search in this area will ascertain the extent to which neutrino transport is instru- 
mental in making a supernova explode. Other bonuses include the elucidation of 
the possible role of supernovae and neutrinos in r— process nucleosynthesis. 

The main issues that emerge from PNS studies concern the metastability and 
subsequent collapse to a black hole of a PNS containing quark matter, or other 
types of matter including hyperons or a Bose condensate, which could be ob- 
servable in the v signal. However, discriminating among various compositions 
may require more than one such observation. This highlights the need for break- 
throughs in lattice simulations of QCD at finite baryon density in order to unam- 
biguously determine the EOS of high density matter. In the meantime, intriguing 
possible extensions of supernova and PNS simulations with npQ and npK matter 
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include the consideration of heterogenoeus structures and quark matter superflu- 

idity (H). 
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Figure 1: Neutrino mean free paths in matter with nucleons only (left panels). 
Right panels show ratios of mean free paths in matter without and with hyperons. 
Top panels show scattering mean free paths (common to all neutrino species). 
Bottom panels show v e mean free paths including absorption reactions. The 
figure is from Ref. (Jl6|) . 




Figure 2: Left: Charged current inverse neutrino mean free paths. Right: 
Comparison of scattering mean free paths in neutrino poor matter at fixed entropy 
in matter containing nucleons and also hyperons. The figure is from Ref. (16). 
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Figure 3: The density and temperature dependences of the neutral current 
mean free paths for neutrino free matter in the field-theoretical model GM3. The 
upper left panel shows the Hartree results for the case E v = 3T. The influence 
of the spin correlations introduced via the Migdal parameter </ is strong, as can 
be deduced from the results shown in the upper right and bottom panels. The 
figure is from Ref. (|l7|). 
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Figure 4: The density and temperature dependences of the charged current neu- 
trino mean free path in (5— stable matter for the GM3 model assuming Yl = 0.4. 
Results for the Hartree approximation (left panel) are compared with those in- 
cluding RPA corrections (right panel) with g' = 0.6. The figure is from Ref. 
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Figure 5: v e mean free paths from various particles in matter containing a mixed 
phase of quarks and hadrons. Thick lines show the extent of the mixed phase 
region. Left (right) panels show scattering (absorption) mean free paths. The 
upper (lower) panels correspond to the neutrino-trapped (neutrino-free) era. Ver- 
tical dashed lines labelled u\a and u max indicate the central densities of 1.4 M Q 
and maximum mass (2.22 M© for the upper panels and 1.89 M for the lower 
panels) stars, respectively. 
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Figure 6: Transport neutrino mean free paths in heterogeneous media. Solid lines 
are for matter in a mixed phase containing kaons (left panel) and quarks (right 
panel), and dashed curves are for uniform matter. This figure is from Ref. (23). 
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Figure 7: Neutrino asymmetry growth curves driven by v T <-> v s and v T «-> v s 
oscillations. The mixing angle is selected to be sin 26> = IO -8 . The three curves 
correspond to Am 2 = —0.5, —50 and —5000 eV 2 , reading from left to right. This 
figure is taken from Ref. (|49|). 
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Figure 8: The main stages of evolution of a neutron star. Shading indicates, 
approximately, relative temperatures. This figure is from Ref. (p5|). 



Neutrino Propagation in Dense Astrophysical Systems 



73 




Figure 9: The evolution of the average energy and total luminosity of neutrinos 
in PNSs composed of baryons only (left panel) and baryons and hyperons (right 
panel). The figure is from Ref. (|8q). 
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Figure 10: The upper panel shows the total emitted neutrino luminosity for the 
evolution of a PNS using the opacities of Ref. (fT^). Here u = ns/n®. The lower 
panel shows the ratio of the luminosities obtained for models with correlation 
corrections to the baseline (Hartree approximation) model. This figure is from 
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Figure 11: Emitted neutrino luminosity in PNSs (n = ns/no). This figure is 
fromRef. 
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Figure 12: The evolution of the total neutrino luminosity for npQ PNSs. Shaded 
bands illustrate the limiting luminosities corresponding to a count rate of 0.2 
Hz, assuming a supernova distance of 50 kpc for 1MB and Kamioka, and 8.5 kpc 
for SNO and SuperK. The widths of the shaded regions represent uncertainties 
in the average neutrino energy from the use of a diffusion scheme for neutrino 
transport. This figure is from Ref. (|90|). 
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Figure 13: Lifetimes of metastable stars versus the PNS baryon mass Mb- Thick 
lines denote cases in which the maximum gravitational masses of cold, catalyzed 
stars are near 1.45 M , which minimizes the metastability lifetimes. The thin 
lines for the npQ and npH cases are for EOSs with larger maximum gravitational 
masses (1.85 and 1.55 M , respectively.) This figure is from Ref. (pQ). 



